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Abstract. We discuss how we can obtain new quantum superintegrable Hamiltonians allowing 
the separation of variables in Cartesian coordinates with higher order integrals of motion 
from ladder operators. We also discuss how higher order supersymmetric quantum mechanics 
can be used to obtain systems with higher order ladder operators and their polynomial 
Heisenberg algebra. We present a new family of superintegrable systems involving the fifth 
Painleve transcendent which possess fourth order ladder operators constructed from second 
order supersymmetric quantum mechanics. We present the polynomial algebra of this family of 
superintegrable systems. 



1. Introduction 

The most well known superintegrable systems are the Kepler-Coulomb system and the harmonic 
oscillator [1]. Superintegrable systems with second order integrals of motion in two-dimensional 
[2] and three dimensional Euclidean spaces [3,4] were classified. However, a systematic search for 
superintegrable systems with third order integrals of motion in classical and quantum mechanics 
is more recent and was started in Ref.5. Classical and quantum systems with a second and a 
third order integrals of motion allowing separation of variables in Cartesian coordinates were 
obtained by S. Gravel [6] and five of the 14 quantum systems were written in terms of the first, 
second and fourth Painleve transcendents [7]. These quantum systems were studied [8-10] from 
the point of view of supersymmetric quantum mechanics and their polynomial algebras. For 
a review of two-dimensional superintegrable systems we refer the reader to the following paper 
[11]. In recent years, many articles were devoted to integrable and superintegrable systems with 
higher order integrals of motion [12-31]. 

A relation between first and second order ladder operators and integrals of motion of 
superintegrable systems were acknowledged by many authors [1-4,10,15,16,32-34]. A relation 
between first order supercharges [35] and integrals of motion were also discussed [36-39]. Recently, 
we pointed out [8-10,12,13] how superintegrability with higher order integrals, higher order 
ladder operators and higher order supersymmetric quantum mechanics [40-44] are related. We 
discussed [10] how from known one-dimensional Hamiltonians that possess ladder operators 
we can form multi-dimensional superintegrable systems, their integrals and their polynomial 
algebra. We constructed a new infinite family of superintegrable systems involving the fourth 
Painleve transcendent from third order ladder operators. We presented a classical analog of 
this construction and obtained an infinite family of classical superintegrable using third order 
ladder operators [12]. The factorization of Hamiltonians in supersymmetric quantum mechanics 



is not necessarily unique and can be used as a tool for obtaining from known systems new 
superintegrable systems with higher order integrals of motion [13]. 

The purpose of this paper is to discuss how from systems allowing two intertwining relations 
with second order supercharges [44], which are related to fourth order ladder operators, we can 
generate an infinite family of superintegrable systems involving the fifth Painleve transcendent. 
In Section 2, we present a method to construct superintegrable systems introduced in Ref.10. 
In Section 3, we discuss the relation between supersymmetric quantum mechanics and ladder 
operators. In Section 4, we recall results of Ref.44 and using these results we present explicitly 
the system with fourth order ladder operators, its ladder operators and zero modes in terms 
of the fifth Painleve transcendent. In Section 5, using results of Section 2 and Section 4, we 
obtain a new family of superintegrable systems and we present their integrals of motion and the 
polynomial algebra generated by these integrals. 

2. Superintegrability and ladder operators 

Let us consider the following two-dimensional Hamiltonian allowing separation of variables in 
Cartesian coordinates: 

H = H 1 + H 2 = ^- + ^- + V 1 (x 1 ) + V 2 (x 2 ). (1) 

We suppose the existence of polynomial ladder operators (oji are constants and Qi polynomials 
of order ki) : 

[H it A\] = m A\, [H, Ai] = -UiAi, (2) 

AiAt = Qi(Hi + u}i), A\Ai = Qi(Hi), i = l,2. (3) 

The polynomial Heisenberg algebra generated by {Hi, Ai, Al} provides information on the 
spectrum of Hi. The annihilation operator can allow at most ki zero modes (i.e. a state such 
Aiip = 0). In each axis we can have at most ki infinite ladders by acting iteratively with the 
creation operator. The creation operator can also allow zero modes and we can have finite 
ladders. 

From these operators, we construct the following integrals of motion (of order 2, k\n\ + k 2 n 2 
and k±ni + k 2 n 2 ) with n\uj\ = n 2 oo 2 = u, n\,n 2 £ N* : 

K = ±-(Ht-H 2 ), H = A^4 n \ I + =A^A™. (4) 

The systems is thus superintegrable. We can also consider the integrals I\ = J_ — I + and 
I 2 = i_ + I + . These results can be used to study known superintegrable systems and also 
construct new ones. Moreover, this construction can be extended for N-dimensional quantum 
systems and also in context of classical systems [12]. These integrals of motion generate a 
polynomial algebra [8-12,45-48] of order (kim + k 2 n 2 — 1) 

[K,I ± ]=±I±, [H,I + ] = F nun2 (K + l,H)-F nun2 (K,H), (5) 

ni tt 712 H 

F nun2 (K,H) = Y[Qi(— + n^K - ( ni - iMjlQ^— - n 2 uj 2 K + ju 2 ). (6) 

i=l j=l 

The equations (5) and (6) form a deformed u(2) algebra and the finite dimensional representation 
modules can be found using a realization as a generalized deformed oscillator algebra [49]. The 
operators = i+, b = I-, N = K — u ( u is a constant) and 3>(iV, u, H) = F ni>n2 (K, H) satisfy 
a deformed oscillator algebra : 



[N,b t ] = b t , [N,b] = -b, b t b = <5>{N,u,H), bb t = $(N+l,u,H) . (7) 



The function <3?(iV) is called the "structure function". Following [7-10,47,48] we have an energy 
dependent Fock space of dimension p+1 if 

$(0,u,£)=0, $(p + l,u,E) =0, $(n,u,E) > 0, n = 1, 2, (8) 

The Fock space is defined by 

H\E,n >= E\E,n >, 2V|E, n >= n\E, n > b\E, 0>=0, (9) 

b t \n>= y/^r+T^E)\E,n + l>, b\n >= \J 3>(n, u, E)\E, n — 1 > . (10) 

The relations given by equation (8) can be used to obtain the finite dimensional unitary 
representations and the corresponding degenerate energy spectrum of superintegrable systems 
[8-11,47,48]. 



2.1. Systems with ladder operators of first, second and third order 

Infinite families of superintegrable systems in two-dimensional Euclidean space from first [1], 
second [15] and third [10] order ladder operators were discussed in previous articles (with 
n\U)\ = H2U2 = io, n\,U2 G N*): 

V(xi,x 2 ) = \^oj\x\ + u) 2 2 x 2 ,), (11) 



V( Xl ,x 2 ) = ^K 2 xi +64^ + 4 + 4' ( 12 ) 

£ X X 2 



V{x l ,x 2 ) = —{xi + \/- ? 4(i)(W^i),ai.Pi) H y~ P WV -fi X i> a i>Pi) + ~vei -«i) 

(13) 

u 2 I / w 2 « n\2 (-2^2 ' 1^2 n ^2, 

+yO r 2 + y— P 4 (2)(W Y^2,a2,P 2 )) + — g— P 4( 2)( Y -7-^2,02, /3 2 ) + — (e 2 -a 2 ) 

= 9 p W ^ + + 4zi? (i) (zi) + 2(z 2 - Oi)P4(i)K) + TT^TTT' ( 14 ) 

-P 4 (i)(^) = P 4(i)(^ 5 «i)ft)) Z% = \f^ Xi - 

The integrals of motion, the polynomial algebras , the finite-dimensional unitary repsentations 
and the corresponding degenerate energy spectrum were given respectively for superintegrable 
systems with first order ladder operators in Ref. 33 and systems with second and third order 
ladder operators in Ref. 10. 



3. Supersymmetry and ladder operators 

Let us recall some results of first and second order super symmetric quantum mechanics (also 
known as intertwining or factorization method) [35] 
Let us consider the following intertwining relation 

H a tf = L^H b , (15) 

with a first order intertwining operator (also called supercharge) and two Hamiltonians of the 
form 

h 2 d 2 Id 
#i = -y^2 +^(x), = —(-h—+a(x)), i = a,b. (16) 



These relations lead to 



V b = V a - ha', 2aV b - h 2 a" = 2aV a - 2W' a . (17) 

We obtain by putting the first relation in the second and integrating the result 

ha' + a 2 = 2(V a (x) - e). (18) 

These equations are equivalent to the factorization method 

H a = L^L + e, H b = W + e. (19) 

If one of these systems is exactly solvable and allows ladder operators, we can obtained the 
wavefunctions (up to a possible zero mode state), the ladder operators and the energy spectrum 
of its superpartners using the supercharges. The operators Aa and A\ (A a and A b ) respectively 
the raising (lowering) operators of the Hamiltonians H a and H b are related by A a = l} A\L 
(A a = L'A b L). Supersymmetric quantum mechanics can thus be used to generate systems with 
higher order ladder operators. 

Let us now consider two Hamiltonians [40-44] 

L] = V h2 ^- h9{x) T x +h{x)) - (20) 

We obtain the following three equations 

V b = V a - hg(x)', (21) 

tf&L _ Hh{x)' - g(x)V b = 2W' a - g(x)V a , (22) 

2h(x)V b - h 2 h{x)" = 2h 2 VQ - 2hg{x)V' a + 2h(x)V a . (23) 
From the equations (21), (22) and (23) we can obtain 

v ^ {x) = h Ur) - h %k ±hA r + A r + Th + 2SWH' (24) 

h(x) = -h 29 ^ + h 2 ^-h^ + 9 -^f- * (25) 
v ' 2g(x) Ag(x) 2 2 4 g{x) 2 h v ' 



Moreover, we have 



The ladder operators of the Hamiltonians H a and H b are thus related by the second order 
supercharges. 

A first and second order supersymmetry can also be combined [41] to obtain systems with 
third order ladder operators related to superintegrable systems obtained by S. Gravel [6,8-10]. If 
we consider a first H a L\ = L\(H b + hw) and a second order intertwining relations HaL^ = L\H b , 
we can construct the following third order supercharges A* = l\l<2 and A = L^Lx- They satisfy 
respectively the following relations : [H a , A'] = hujA^ and [iJ a ,A] = —hwA. These operators 
can thus be interpreted as raising and lowering operators. These potentials can be written in 
terms of the fourth Painleve transcendent. The equation (13) consists in a sum of two such 
one-dimensional systems. 



4. Systems with fourth order ladder operators 

Let us recall results from Ref.44 and consider two Hamiltonians 

H i = ~Y^ + Vi ^ i = aA (27) 

intertwined by two second order supercharges L\ and L 2 

L{ = \{f?^- 2 - h 9l {x)j- x + hi(x)), 4 = i(n 2 ^ - hg 2 {x)j- + h 2 {x)). (28) 
The interwining equations have the following form 

H h L\ = L\H a , (29) 

H a L{ = L\{H b + hw). (30) 

We construct from these supercharges the following fourth order ladder operators of the 
Hamiltonian H a and their polynomial Heisenberg algebra 

^ = l\l\, A = L 2 L x , (31) 

[H ai = huA\ [H a , A] = -HujA, (32) 
4 

A* A = Q(H a ) = Y\{H a - ej), AA^ = Q(H a + hw), (33) 
3 

with ei j2 = huj + || ± y^H") ^3,4 = fl =•= The operators (A^ = L\l\ and A = L\L 2 ) are 

ladder operators of the Hamiltonian H 2 . 

The intertwining relations (29) and (30) give six equations. Using the following changes of 
variables z = y/j-x, gi = \ffajjg~i, hi = uihhi (with i=l,2), Vk = fuoVk (with k=a,b) and e'j = ^ 
(with j = 1, 2, 3, 4) we obtain 

r7 " /r 1 n 2 (J J \i 

91 ( 91 s2 _ ~ / , 91 , l e 3 e A) 1 J 1 J _ 9T/ ~ (o A \ 

TT^r - y-z-r) T Si + — : — I h e 3 + e 4 - 2K a6 , (34) 

2^1 2gi 4 3^ 



n " n ' ~2 (J J \2 
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(35) 



f^ + ^-^ + ^-^^lP^^ (36) 
52" 

V6 = y a + 5 Y-l = K-<72'. (38) 

From the equations (34), (35) and (38), we can obtain the following nonlinear second order 
differential equation [44] 

Sk" = 9 ( ~ 2g l + ^\ (gV) 2 + -, 9 \ , 9~2' + +[2zg~ 2 {-g~ 2 + z)r 1 [-2zg~ 2 5 + {hz 2 + ^ + A)g^ (39) 
2g 2 {-g 2 + z) z{-g 2 + z) 

-iz(z 2 + 4e + 2)5 2 3 + [z A + 4(2e + l)z 2 + 4(A 2 2 - A?) - l]^ 2 - 4A 2 ^(-2g 2 + z)} 

with ( e = ±(e' 3 + e' 4 ) - ^ + e' 2 ), A x = e' x - e' 2 , A 2 = e 3 - e 4 ). 

The equation (39) can be transformed into the fifth Painleve transcendent by considering 

{92 = w=i and y = z2 ) 

d 2 VY_ ( _}_ J- )( dW s.2 ldW {W -I) 2 aW 2 + b cW dW{W + 1) 

A 2 A 2 1 1 

with the following parameters a = - 2 X , b = — ^, c = — e — 5, and d = — g. 



4-1- Systems in terms of the fifth Painleve transcendent 

Using these results, we present the Hamiltonians H a and Ht>, the supercharges Li, L2, the ladder 
operators L and and their zero modes in terms of the fifth Painleve transcendent. 
The potential V(x) = V a (x) — \ {e\ + £2) is given by the following expression 

_ _ * . w 2 4(P 5 + P') 2 - P 2 2 /j 2 1 b-aPh t ^ (l + 2cP 5 ), 

= y (1 + > 5 -im )x + " 6 " s + " M1 + WiF' (41) 

with (P 5 = p 5 (^x 2 ,a,6,c,-|)). 

The functions hi,h2,gi and 32 in terms of the fifth Painleve transcendent are given by 

[U 2(4c-P 5 ) N fc2 , -8fe + P 5 (l-8a + 8fe + 8aP 5 ) , (A0 , 

9i = p^i + n x ' /ti=M 4 ( -i+p 5 ) )+M } (42) 

^ 2, x 2 (P 5 (-l + (-1 + P 5 )P 5 ) + 4P 5 P^ - 4P^ 2 ) 
+W ( 4(-l + P 5 )2F 5 )' 

^ = 2(1 + 4cP 5 ) -86 + P 5 (l - 8a + 86 + 8aP 5 ) 

y P 5 -l V 4(-l + P 5 ) ; V 4x 2 P 5 

2 -x 2 (P 2 + P| + 4P* + P 5 (-l + 4PQ) 
+U ( 4(-l + P 5 ) 2 P 5 } - 

There are four possible zero modes for the annihilation operator 

4% oc (±^(-1 + P 5 ) - (1 + c) + ri 4 ( li + p 5 ) ) g/ ( 44 ) 

V> 3) 4 0ce Jl 2*(-i+p 5) p 5 + 2.^ (45) 

with the corresponding energy levels 

E^ 12 = E i%i = hw{-c - 1 T y^)i (46) 

and also four possible zero modes for the creation operator 

-i-^x 2 ±v^:(-i+p5) 2 +p 5 -2^ : i ; |p i ; 
0i ;2 oc e J M-i+flO , (47) 

03,4 cx ( T y -(-1 + P 5 ) + + c) - A I pT^rf Z ) e ^Fi+fe (48) 

with the following energy levels 

'a 
2 



% s =±fc(, --l), ^ 3i4= ^;(_ c __ TA /_). (49) 




Thus, the system can have one, two, three or four infinite sequences of levels and finite sequence 
of levels (singlet, doublet or triplet states) depending of the parameters (a,b,c and oS). A more 
detailed study of these cases should be performed and we let it to a future article. The fifth 
Painleve transcendent also allows rational solutions [50] depending of the parameters. 



5. Superintegrable systems 

Using results from Section 2, we can construct an infinite family of two-dimensional 
superintegrable systems (with n\U)\ = ri2UJ2 — A, ni,ri2 € N*) 

Y 2 8 (p m -iyp 5[l) ,l xf l 1 s p m > A 2{p m -i) h 

(50) 

with {P$(i\ = -P5(j)(^x?, a,i, bi, Ci, — |) i=l,2). The integrals of motion are given by (using (28), 
(31), (42), (43)) 

K=±(H 1 -H 2 ), I_ = A n s4 n \ I + = A\ n 'Af. (51) 

These integrals of motion are of order 2 and 4(/si + fo)- The structure function of the polynomial 
algebra (of order 4(fci + ^2) — 1) is given by 

ni H 712 H 

$>(N, u, H) = Yl Qi(— + n^N + u) - (ni - i)wi) Q 2 {— - n 2 0J 2 {N + u) + juj 2 ), (52) 

i=i j=i 



Q i m = l6l[(H l -E ( £), 4) 2 = T ^ n /|, E f 3A = fiuHi-ct - 5 T (53) 
fe=i 



6. Conclusion 

We presented a new family of superintegrable systems with higher integrals of motion from 
results obtained in context of supersymmetric quantum mechanics [44] using a method discussed 
Ref.10. The integrals of motion were constructed from the fourth order ladder operators that are 
product of second order supercharges. We presented also the polynomial algebra generated by 
these integrals of motion. The classical case of systems allowing fourth order ladder operators 
would differ of the quantum case. We let to a future article to obtain and study such systems 
and the corresponding classical superintegrable system. 

The relation between superintegrability, ladder operators and supersymmetric quantum 
mechanics can be used to study known superintegrable systems or obtain new superintegrable 
systems with higher order integrals of motion [10,13,39] and need to be more explored. 

The general equations for third order interwining relations for stationary Schrodinger 
operators [51] and second order intertwining relations for non-stationary Schrodinger operators 
[52] were obtained and may be used to construct superintegrable systems. 

Moreover, systems written in terms of the Painleve transcendents appear in many contexts in 
quantum mechanics [6-11,13,14,52-58] and recently a superintegrable system involving the sixth 
Painleve transcendent was obtained [14]. The results obtained in this paper also extend the 
number of known superintegrable systems involving the Painleve transcendents. 
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